Transverse-field Ising spin chain with inhomogeneous disorder 
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We consider the critical and off-critical properties at the boundary of the random transverse- 
field Ising spin chain when the distribution of the couplings and/or transverse fields, at a distance 
I from the surface, deviates from its uniform bulk value by terms of order l~ K with an amplitude 
A. Exact results are obtained using a correspondence between the surface magnetization of the 
model and the surviving probability of a random walk with time-dependent absorbing boundary 
conditions. For slow enough decay, k < 1/2, the inhomogeneity is relevant: Either the surface 
stays ordered at the bulk critical point or the average surface magnetization displays an essential 
singularity, depending on the sign of A. In the marginal situation, k = 1/2, the average surface 
magnetization decays as a power law with a continuously varying, A-dependent, critical exponent 
which is obtained analytically. The behavior of the critical and off-critical autocorrelation functions 
as well as the scaling form of the probability distributions for the surface magnetization and the 
first gaps are determined through a phenomenological scaling theory. In the Griffiths phase, the 
properties of the Griffiths-McCoy singularities are not affected by the inhomogeneity. The various 
results are checked using numerical methods based on a mapping to free fermions. 

05.50.+q, 64.60.Fr, 68.35.Rh 



I. INTRODUCTION AND SUMMARY 

Quantum systems with quenched disorder have re- 
ceived much attention recently due to their unusual 
static and dynamical properties.™ Many of these features 
can be observed in one-dimensional models for which sev- 
eral exact and conjectured results are available in the case 
of homogeneous disorder. Among these models we shall 
consider the random trans verse- field Ising model (RTIM) 
defined by the Hamiltonian 



The homogeneously disordered model, with position- 
independent distributions tti(J) = 7r(J) and pi(h) = 
p(h), has received much attention, especially after Fisher 
has obtained new striking results about static critical 
properties and equal-time correlations, using a renorma 
ization group method based on a decimation procedure. 
Many of Fisher's results have been later verified by nu- 
merical methods and some other analytical, scaling and 
numerical results have been obtained concerning the dy- 
namical properties of the RTIM-,a<id the behavior of dif- 
ferent probability distributions u 



(i.i) 



For the homogeneous model in (1.1) one defines the 
quantum control parameter as 



where af and af are Pauli matrices at site I. The ex- 
change couplings Ji and the transverse fields hi are ran- 
dom variables with probability distributions 7T;(J) and 
pi(h), respectively, which are independent from each 
other but may depend on the position Note that in 
one dimension the couplings and the fields can be taken 
positive without loss of generality. 



6 = [ln/i]av - [In J] a 



(1.2) 



where [. . .] av denotes an average over quenched disorder. 
The system is in the ferromagnetic (paramagnetic) phase 
when the couplings are in average stronger (weaker) than 
the transverse fields, thus for 5 < (5 > 0). The critical 
point corresponds to 6 = 0. 

In this paper we concentrate on the surface proper- 
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ties of the RTIM. The average surface magnetization, 
[m s (5)]avj which can be studied using a mapping to a 
problem of surviving random walks (RW),H behaves as 
[m s (£)] av ~ —8 close to the critical point, whereas the 
correlation length diverges as [£] av ~ |<5| 2 . Thus the cor- 
responding critical exponents for the average quantities 
are: 

13. = 1, v = 2. (1.3) 

The scaling is strongly anisotropic at the critical point of 
the homogeneous RTIM: The time scale set by the relax- 
ation time r r and the length scale £ are related through 

lnr r ~^/2. (1 . 4) 

As a consequence, the imaginary-time spin-spin autocor- 
relations are logarithmically slow at the critical point 
with 

[G s (r)] av ~ (lnr)" 1 , 8 = 0, (1.5) 

for the surface spins. In the Griffiths phase, on the para- 
magnetic side of the critical point, the autocorrelations 
are still anomalous and decay like a power, 

[G(r)]av ~ t- x '*W , 0<8<8 G , (1.6) 

with a dynamical exponent zJJi) given by the positive 
root of the following equations 




As shown recently,Q all the singular quantities in the 
Griffiths phase (susceptibility, non-linear susceptibility, 
energy-density autocorrelations, etc) involve the dynam- 
ical exponent z{8). 

In many physical systems the disorder is not homoge- 
neous. For example a free surface or an internal defect 
line may locally modify the distribution of the couplings 
and/or fields. Here we consider surface induced inhomo- 
geneities which are characterized by a power-law varia- 
tion in the probability distribution: ni{J) — tt(J) ~ l~ K 
and/or pi(h) — p(h) ~ l~ K , for I 3> 1, such that the local 
control parameter in ( |l.2| ) varies as 

8(l) = 8-Al~ K . (1.8) 

Our choice for the functional form of the inhomogeneous 
disorder is analogous to the variation of the couplings in 
the so-called extended surface defect problem first in- 
troduced by Hilhorst and van Lceuwen (HvL) for the 
two-dimensional classical Ising model.u This type of in- 
homogeneity has been later studied for other models and 
different geometries. For a review see Ref. 

Informations about the critical properties of the RTIM 
with inhomogeneous surface disorder have been obtained 
by combining different approaches: The calculation of the 
surface magnetization can be mapped onto a problem of 
surviving RW's, allowing us to deduce some exact results. 
The asymptotic behavior of the surface autocorrelation 
function and the form of the probability distribution for 
different quantities follow from a phenomenological scal- 
ing theory. They have been checked through large scale 
numerical calculations using the free fermion formulation 
of the RTIM. 



TABLE I. Summary of the surface critical properties of the inhomogeneous RTIM. 
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-/»«(-*) 



1-2k 

exp[— const (In r) 1 ~ K 



Our main results are summarized in Table I. The sur- 
face critical behavior of the inhomogeneous RTIM de- 
pends on the value of the decay exponent k. For fast 
enough decay of the inhomogeneity, i.e., for n > 1/2, 
the surface critical properties of the model are the same 
as for the homogeneous RTIM. Specifical ly, t he critical 
expon ents keep the values given in Eq. (1.3) and rela- 
tions (|l.4[ ) and (1.5) remain valid. 



In the borderline 1/2, the effect of the in- 

homogeneity is marginal: the average surface magne- 
tization has a distribution-dependent critical exponent, 
/3 a (A), which has been calculated analytically. This ex- 
ponent governs also the logarithmic decay of the sur face 
spin-spin autocorrelation function, like in Eq. (1.5) for 
the homogeneous case, but with a different value. The 
scaling relation between r r and £ keeps the form given in 
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Eq. ( |l.4| ) for the homogeneous model. The same is true 
of the scaling behavior of the probability distributions 
for the energy gap e, and the surface magnetization m s , 
at large system size L. The scaling function of the lat- 
ter, however, has a different asymptotic behavior, which 
is connected to the different (^-dependence of the average 
surface magnetization in the two cases. 

For a slowe r de cay, such that k < 1/2, the scaling re- 
lation in Eq. fll.4| ) is modified and leads to a new scaling 
form for the probability distributions of the energy gap 
and the suface magnetization. For enhanced surface cou- 
plings, A > 0, the average surface magnetization remains 
finite at the bulk critical point. It vanishes as a power 
of A when A goes to zero from above. For weakened 
surface couplings, the average surface magnetization dis- 
plays an essential singularity in 8 instead of a power law, 
whereas the surface autocorrelation function has an en- 
hanced power-law decay. 

In the numerical calculations we used two types of dis- 
tributions for the disorder. In the inhomogeneous binary 
distribution, the couplings take the values A > 1 and 
A" 1 with probability pi = |(1 + Abl~ K ) and qt = 1 — pi, 
respectively, while the transverse field remains constant: 



P {h) 



■ piS(J - A) + qiS(J — A ) 
S(h-h ). 



(1.9) 



According to Eq. (1.2) ho = 1 at the b ulk critical point 
and the local control parameter in Eq. (1.8) involves the 
parameter A = A^lnA. In the bulk Griffiths phase, 
1 < hp < A, the dynamical exponent, which follows from 
Eq. (1.7), is the solution of the implicit equation: 



cosh 



In A 



(1.10) 



In the uniform distribution, both the couplings and the 
fields have rectangular distributions: 



n(J) = ir(J) = 

o,m - { if"'" 1 " 



1 if < J < 1 

otherwise 

1 ifO<fc<fto(Z) 
otherwise 



(1.11) 



where the inhomogeneity now affects the distribution of 
the fields with ho (I) = hp — A u l~ K . The bulk critical point 
is still at ho = 1 whereas A — A u in the expression of 
the local control parameter in Eq. (1.8). The dynamical 
exponent follows from 



:ln(l - 



= - In h . 



(1.12) 



and the domain of the Griffiths phase now extends to 
1 < ho < oo. 

The structure of the paper is the following: In Sec. II 
we use the free fermion formulation of the RTIM to ex- 
press different surface quantities. The relation between 
the surface magnetization of the inhomogeneous RTIM 



and an absorbing RW problem is treated in Sec. III. It is 
exploited there to develop a relevance-irrelevance crite- 
rion. Our results are presented in Sec. IV and discussed 
in Sec. V. 



II. FREE-FERMIONIC EXPRESSION OF 
SURFACE QUANTITIES 

We start by considering the imaginary-time surface au- 
tocorrelation function of the RTIM 



G s (r) = (0K(tK(Q)|0) 

KzKlO^exphr^-^o)] 



(2.1) 



where |0) and \i) denote the ground state and the i-th 
excited state of Tt in Eq. (1.1) with energies Eq and Ei, 
respectively. With symmetry-breaking boundary condi- 
tions, i.e., with a fixed spin at the right end of the chain, 
cr£ = ±1, the ground state of the system is degenerate 
and the autocorrelation function asymptotically behaves 



as Iim T _ ) .oo G s (t) 
tion is given by: 



so that the surface magnetiza- 



(iK|o) 



(2.2) 



In order to calculate m s and other quantities we use 
the method of Lieb, Schultz and MattisEil and transform 
TL into a free fermion Hamiltonian, 



9=1 



VqVq 



(2.3) 



where the 77 J 's (?7 g 's) are fermion creation (annihilation) 
operators. For free boundary conditions the fermion ex- 
citation energies e q are obtained through the diagonal- 
ization of the following 2L x 2L tridiagonal matrix:E3 



T = 



( hi 
hi Ji 
Ji h 2 

h 2 



V 



'■■ Jl-i 
Jl-i h L 
h L / 



(2.4) 



The components of the eigenvectors V 9 are written as 
V q {2l-l) = -<p q (l) and V q {2l) = I = 1, 2, . . . , L. 

Changing 4> q into — <p q one obtains the eigenvector corre- 
sponding to — e q , thus we confine ourselves to that part 
of the spectrum corresponding to e q > which contains 
all the needed informationJlj Fixing the surface spin at 
I = L amounts to take h^ = 0. Then ei = and the 
ground state of the Hamiltonian is degenerate. 

The surface magnetization and the spin-spin autocor- 
relation function in the free-fermion description are cal- 
culated using Wick's theorem and can be expressed in 
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FIG. 1. Random walk correspondance of the distribution 
of couplings. 

The distribution of the couplings in a sample can be 
put in correspondence with a RW, y(t), the time axis t 
corresponding to the position j along the spin chain. As 
shown in Fig. [|, the walker starts at y = 0, its j-th step 
is in the positive (negative) y-direction when Jj = A 
(A -1 ) and absorbing boundary conditions are assumed 
for y < 0. As a consequence, the surface magnetization 
of a sample is nonvanishing only when the corresponding 
walk survives until t = L, i.e., if it never crosses the time 
axis during L — 1 steps. Thus the fraction of L-step sur- 
viving walks, Psurv(t)\ t =L ~ L^ 1 / 2 , gives the fraction of 
rare events for which a sample has a nonvanishing sur- 
face magnetization. Since the average value of the surface 
magnetization is determined by the rare events, we have 
the correspondence: [m s (L)] av ~ P S urv(t)\ t=L ~ L^ 1 / 2 
from which the value of the critical e xpon ent x s ~ (3 s /i> = 
1/2 follows, in agreement with Eq. ( |l.3| ). 

Let us consid er th e average of the logarithm of the 
products in Eq. (2.5) which may be written as 



□ 



1 'h.- 2 



3 = 1 



2^([ln/i J -] av -[lnJ,-] av ) = 25r 



(3.1) 



□ 



D 



□ 



using the definition of the control parameter in Eq. fll.2j ). 
In the RW picture this quantity is proportional to the av- 
erage position of the walker after I steps for a free walk 
and vanishes at criticality. In the off-critical situation, 
<5 7^ 0, the average position of the walker grows lin- 
early with which is equivalent to having a nonzero bias 
Sw = Qw — Piv 7^ where p w (q w ) denotes the probability 
that the walker makes a step in the positive (negative) y- 
direction. Thus the correspondence between RTIM and 
RW still applies in the off-critical situation with: 



,(S,L)] e 



r(S w ,t)\ 



t=L 



(3.2) 



From the known behavior of the surviving probability, 
one deduces the values of the RTIM critical exponents 
given in Eq. (1.3). 



4 



For the inhomogeneous RTIM the local quantum con- 
trol parameter has a smooth position dependence which, 
at the bulk critical point, is given by 6(1) = —Al~ K ac- 
cording to Eq. The corresponding RW has a lo- 
cally varying bias with the same type of asymptotic de- 
pendence, 6 w (l) = —A w l~ K . Consequently the average 
motion of the walker is parabolic: 



A,, 



(3.3) 



1=1 



Under the change of variable y(t) — > y(t) — y p (t), the sur- 
viving probability of the inhomogeneously drifted walker 
is also the surviving probability of an unbiased walker, 
however with a time-dependent absorbing boundary con- 
dition at y(t) < —y p (t). 

The surviving probability of a random walker with 
time-dependent absorbing boundaries has already been 
studied in the mathematically and physicalllSO litera- 
ture. In a continuum approximation, it follows from 
the solution of the diffusion equation with appropriate 
boundary conditions, 



d 



d 



-P(y,t) = Dj-^P(y,t) , P[-y p (t),t]=0. (3.4) 



dt 



dy 



Here P(y, t) is the probability density for the position of 
the walker at time t so that the surviving probability is 
given by 



P a „rv(i)= / P(y,t)dy. 
'-»*(*) 



(3.5) 



The behavior of the surviving probability depends on the 
value of the decay exponent re. For n > 1/2, the drift of 
the absorbing boundary in Eq. ( |3.3| ) is slower then the 
diffusive motion of the walker, typically given by 

(Dtfl\ 



Vd(t) 



(3.6) 



thus the surviving probability behaves as in the static 
case. When re < 1/2, the drift of the absorbing bound- 
ary is faster than the diffusive motion of the walker and 
leads to a new behavior for the surviving probability. 
For A w > 0, since the distance to the moving boundary 
grows in time, the surviving probability approaches a fi- 
nite limit. On the contrary, for A w < 0, the boundary 
moves towards the walker and the surviving probability 
decreases with a fast, stretched-exponential dependence 
on t. Finally, in the borderline case re = 1/2 where the 
drift of the boundary and the diffusive motion have the 
same dependence on t, like in the static case the surviv- 
ing probability decays as a power, P surv (6 W , t) ~ t~ e ( A ™\ 
however with a continuously varying critical exponent. 

Now turning back to the problem of the inhomoge- 
neous RTIM, the RW analogy can be used to formulate 
a relevance-irrelevance criterion: According to the pre- 
ceding discussion, the inhomogenous surface disorder is 
a relevant (irrelevant) perturbation for the surface criti- 
cal properties of the RTIM when re < 1/2 (> 1/2) and a 
marginal one when re = 1/2. 



IV. RESULTS 



In this Section we study the surface critical proper- 
ties of the inhomogeneous RTIM, considering the size- 
dependence at criticality as well as the 5-dependence for 
various quantities. We treat successively relevant and 
marginal perturbations and close with a discussion of the 
behavior in the Griffiths phase. 

In the numerical calculations we used up to 5 10 4 real- 
izations of disorder with chain sizes up to L = 2 7 except 
for the surface magnetization for which we used up to 
10 6 realizations with system sizes up to L = 2 14 . 



A. Critical behavior with relevant inhomogeneity 



Let us first consider the probability distribution of 
lnm s on finite samples with length L at the bulk crit- 
ical point 5 = 0. According to Eq. ( |2.8| ), [lnrn s ] av is 
expected to scale as [lnPJ^_ 1 ( Ji/hi)] av when the pertur- 
bation tend s to reduce the surface order, i.e., when A < 
in Eq. (1.8). Thus one obtains: 



[hlTOsJav ~ ^([ln J;] av - [ln/lijav 
(=1 

L 

= A s ^r K ~AL 1 - K . (4.1) 



i=i 



The typical magnetization, define through [In m s 



ln[r 



tsjtypj 



has a stretched-exponential behavior, 



[TO s jtyp ~ exp(const AL 1 " 1 ^ 1 ), when A < 0. For the 
probability distribution of In m s , Eq. j4.l| ) suggests the 
following scaling form: 



P m (In m s ,L) 



-Pn 



In) 



(4.2) 



Indeed, as shown in Fig . the numerical results are in 
agreement with Eq. (|h2|). 

The asymptotic behavior of the scaling function P m (v) 
as v — > depends on the signe of A. For enhanced sur- 
face couplings, A > 0, there is a nonvanishing surface 
magnetization at the b ulk critical point as L — > oo, thus 
the powers of L in Eq. (4.2) must cancel as seen in Fig. || 
and \im v ^Q P m (v) ~ v~ x so that P m (lnm s ) ~ l/lnm s . 
For reduced surface couplings, A < 0, lim^o P m (v) = 0, 
indicating a vanishing surface magnetization. 
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FIG. 2. Probability distribution of lnm s obtained for the 
binary distribution with k = 1/4 for A — 0.4 and scaling plot 
for A = -0.4 (inset) for chain sizes L = 2 10 ,2 U , 2 12 . 

Next we calculate the average behavior of the surface 
magnetization, which is determined by the rare events 
with m s = 0(1). Here we use the RW description of 
Section III. Starting with A ~ A w > 0, for A w <C 1 one 
defines the time scal e t* for which the parabolic and diffu- 
sive lengths in Eqs. (3.3) and ([si]) are of the same order, 

Vp{t*) ~ Vd(t*), such that t* ~ A„ 2/(1 ~ 2k) . The surviv- 
ing probability can be estimated by noticing that, if the 
walker is not absorbed up to t* , it will later survive with 
a finite probability. Thus P SUIV {A W > 0) - (i*)" 1 / 2 ~ 

Aw^~ 2k ^ and the average surface magnetization has the 
same behavior at the bulk critical point: 



\m s 



1 



< A< 1. 



(4.3) 



This result is in agreement with the numerical calcula- 
tions as shown in Fig. [j| The error bars give the precision 
of the extrapolation procedure. 

In the case of a shrinking interval between the walker 
and the absorbing wall, A w < 0, the leading behavior of 
the surviving probability can be estimated by looking for 
the fraction of walks with y(t) > —y p {t) which is given 
by P S urv(t) ~ cxp[— const A^ i 1-2 "]. Thus, for reduced 
surface couplings, the average surface magnetization has 
the following finite-size behavior at the critical point: 



[m s ] av ~ exp (-const A 2 L 1_2k ) , k < -, A < . (4.4) 

Numerical results are shown in Fig. ||. One can notice 
that even for such large system size as L — 2 9 the linear 
asymptotic regime is not completely reached. One may 
notice the different size dependence of the typical an d av - 
erage sur face magnetizations at criticality in Eqs. (4.2) 
and (4.4), respectively. 
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FIG. 4. ln[m s ]av as a function of L 1 ~ 2k for the binary dis- 
tribution with ft = 1/4 and A = —0.4. Th e asymptotic linear 
dependence is in agreement with Eq.(4.4). 
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FIG. 3. Average surface magnetization as a function of 
the inhomogeneity parameter for the binary distribution with 
k = 1/4, extrapolated at infinite size. The dott ed line corre- 
sponds to the asymptotic expression in Eq.(4.3). 



To obtain the (5-dependence of the average surface mag- 
netization in the thermodynamic limit, we first determine 
the typical size of the surface region l s which is affected 
by the inhomogeneity for 5 < and A < 0. It is given by 
the condition that quantum fluctuation (~ 8) and inho- 
mogeneity (Al~ K ) contributions to the energy are of the 
same order, from which the relation l s ~ |A| 1 / K |(5|~ 1 / K 
follows. Inserting L ~ l s in Eq. (4.4), we obtain 



[m s ] av ~ exp[-const \A\^ K \8\-^- 2k ^% k < \,A<0, 

(4.5) 

as quoted in Table I. The numerical verification of this 
result is difficult since, close to the critical point, one has 



G 



to use large finite samples, L ^> l s , for which [m s (S, A)] 
is quite small. 

Let us now study the behavior of the energy gap 
ei(L) ~ r" 1 . Here we make use of the relation (2. 



between the gap and the surface magnetizations, from 
which the typical finite-size-scaling behavior lnei(L) ~ 
L 1 ~ K follows, as for the surface magnetization. Thus the 
relation between the time and length scales is 



and the probability distribution of the gap is 



P £ (lnei,L) 



1 



L 1 



-P r 



lnei 



(4.6) 



(4.7) 



to be compared to the distribution of the surface magne- 
tization in Eq. ( |4.2| ). Numerical results shown in Fig. |^ 
are in agreement with this scaling result. 
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FIG. 5. Scaling plot of the probability distribution of the 
gap for the binary distribution with A = 0.5 and k = 1/4. 

In order to obtain an estimate for the average surface 
autocorrelation function one may notice that the disor- 
der being strictly correlated along the time axis, a sample 
with a finite surface magnetization m s = O(l) has also 
a nonvanishing surface autocorrelation function G s (r) ~ 
rrig. Since the fraction of rare events are the same for the 
two quantities, the scaling behavior of [G s (L,r)] av can 
be deduced from the corresponding relations for the av- 
erage surface magnetization. F or e nhanced surface cou- 
plings, A > 0, according to Eq. (pf , lim r ^ 00 [G s (r)] av ~ 
[m s (A)] av — yl 1 /( 1 - 2K ) a t criticality. For reduced sur- 
face couplings, A < 0, the finite-size critical behav- 
ior follows from Eq. (4.4) with lim T ^ 00 [G s (i, r)] av ~ 
[m s (L,A)] av ~ exp[— const A 2 L 1_2k ]. Now, using the 
scaling relation (POI), we obtain 



[G s (r)} s 



exp 



-const A 2 (lnr)" 



(4.8) 



in the thermodynamic limit. This enhanced power-law 
decay is consistent with the numerical results, as shown 
in Fig. |. 
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FIG. 6. Average surface autocorrelation function for the 
uniform distribution with k = 1/4 and A — —0.5. The solid 
straight line is a guide for the eye. 



B. Critical behavior with marginal inhomogeneity 



In the borderline CctS6, — 1/2, several results in 
the previous Section have to be modified. We start 
with the analysis of the average surface magnetization. 
For the extreme binary distribution, with A — * oo in 
Eq. ([0|), there is a one-to-one correspondence between 
samples with a nonvanishing surface magnetization and 
RW's with surviving character, as explained in Section 
III. Thus for this distribution one can deduce f3 s from 
the corresponding behavior of the surviving probability 
in the RW problem. The differenti al eq uation with atu 
sorbing boundary conditions in Eq. (3.4) can be solvcdliZl 
for n — 1/2 in_terms of parabolic cylinder functions of 
order is, D^(x)£3 For zero bulk bias, S w — 0, the surviv- 
ing probability has the asymptotic dependence P SU rv{t) ~ 
t~ e and the exponent 8 is such that 



D 2e (-2A V 



0. 



(4.9) 



taking D = 1/2 for the diffusion constant in Eq. (p.4| 
In the limiting cases the solution takes the form 



A;. 



(4.10) 



for \A W \ -C 1 while for large values of A w one obtains 
asymptotically : 



A w exp[-2Al] 



(4.11) 



Then the correspondence of Eq. (3.2) leads to the scaling 
dimension of the average surface magnetization, 
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= ^ = 0(4 

V 



(4.12) 



where the average correlation length exponent 1/ = 2 and 
= 4 for the extreme binary distribution. 
For other i nhom ogeneous distributions, like in 



Eqs. (1.9) and (1.11), one has to find the relation be- 
tween the parameter of the walk A w and that of the 
RTIM A. The argument goes as follows: At the bulk crit- 
ical point, we compare the integrated change in the con- 
trol parameter induced by the perturbation at a length 
scale I, A(l) = 5Z,- =1 [hi Jj — lnhj] av to its fluctuation 
characterized by T 2 (l) = £^ =1 [(ln Jj — ln/ij) 2 ] av . The 

dimensionless ratio A(l)/T(l) is equal to 24 and y/2A u 
for the binary and the uniform distribution, respectively. 
Since it does not depend on the value of A for the binary 
distribution and since A w — 4 when A — > oo, we can 
identify this ratio to 2A W which gives: 



4;, — 
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1.0 



0.0 



4 binary distribution 

4/v2 uniform distribution 



• exact 
□ binary 
A uniform 



(4.13) 
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FIG. 7. Surface magnetic exponent as a function of the 
inhomogeneity parameter A w . The line is a guide for the eye. 

The surface magnetization exponent thus depends on 
the distribution, however only through the value of the 
parameter A w entering in Eq. (4.12), the functional form 



remaining the same. The average magnetization has been 
obtained numerically for the binary and the uniform dis- 
tribution using 2 10 5 realizations of disorder for chain 
sizes L = 2 5 to 2 14 . The exponents x s m were deduced 
from an extrapolation of two-points approximants using 
the BST algorithms. The results are shown in Fig. [?]. 
The exa ct points in the figure are obtained by solving 
Eq. (4.9) for integer and half-integer values of 9 where 
the parabolic, cylinder functions are related to Hermite 
polynomialsE3. The corresponding A w values are then 
given by the highest negative zeros of the Hermite poly- 
nomial, He2#(— 24u) = 0. For x s m = 0, A w is shifted to 



+oo as it can be seen in Eq. (4.11). Thus the system 
does not display surface order at criticality as it is the 
case for the HvL model for sufficiently enhanced surface 
couplings. 

Next we discuss the properties of the probability dis- 
tribution of In m s . Repeating the argument used above 
for rele vant perturbations, we arrive to the scaling form 
of Eq. (12), however with k — 1/2. We note that the 
same scaling form remains valid for irrelevant pertur- 
bations, with k > 1/2, but the scaling function P m (v) 
has different limiting behaviors when v — > in the two 
cases. While for n > 1/2 it approaches a constant, in 
the marginal case it goes to zero as Unit, _> o P m (v) ~ 
v~ 2x ™( A } +1 . In this way one obtains the proper scal- 
ing behavior for the average surface magnetization. The 
numerical results in Fig. are in agreement with this 
limiting behavior. 



&h 1.0 




L m \ln m I 



FIG. 8. Scaling plot of the probability distribution of In rn s 
for the binary distribution with k = 1/2 for chain sizes 
L = 2 10 ,2 n ,2 12 . 

For the distribution of the energy gap at the criti- 
cal point, one can also repeat the reasoning used for 
the relevant case and both the distribution function in 



Eq. (4.7) and the scaling relation in Eq. (4.6) stay valid, 
with k = 1/2. 

The asymptotic behavior of the average surface auto- 
correlation function can be determined at criticality by 
scaling considerations like in a previous work for the ho- 
mogeneous cascfl Here we make use of the fact that, as 
already explained for relevant perturbations, the aver- 
age autocorrelation function has the same scaling prop- 
erties as the average surface magnetization. Under a scal- 
ing transformation when lengths are rescaled by a factor 
b > 1, such that I' = l/b, the average surface autocorre- 
lation behaves as 



[G.(lnr)U = 

where we used Eq. 
scales. Taking now b 



>[G s (lnr/& 1 / 2 )] 8 



(4.14) 



J) to relate the time and length 
(lnr) 2 we obtain: 



S 



[G s (t)U ~ (lnr) 



-2.r" 



(lnr) 



-0s(A) 



(4.15) 



The last e xpression follows from the exponent relation in 
Eq. (4.12) where the average correlation length exponent 
is v = 2. This scaling behavior has also been verified by 
numerical calculations as shown in Fig. [)| 
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FIG. 9. Average surface autocorrelation function for the bi- 
nary distribution in the marginal case with A = 0.5. Th e soli d 
straight line gives the expected scaling behavior in Eq. (4.15). 



distribution of the second gap is also insensitive to the 
presence of the inhomogeneity and the value of the expo- 
nent z'(5), deduced from the slope, is in good agreement 
with the analytical result in Eq. (1.7) with A = 0. 
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FIG. 10. Probability distribution of the first and second 
gaps in the Griffiths phase with the uniform distribution at 
ho — 3 and A = 0.5 for different values of k. The asymptotic 
slopes of the curves, obtained by least square fits, are given 
in parentheses. 



C. Griffiths phase 

In the Griffiths phase, on the paramagnetic side of the 
critical point, the dynamical properties of random quan- 
tum systems are anomalous. As a consequence of the so- 
called Griffiths-McCoy singularities^ the different dy- 
namical properties can be characterized by one single 
parameter, the dynamical exponent z(5), which is a con- 
tinuous function of the quantum control parameter given 
in Eq. (1.7) for the homogeneous RTIM. Here we calcu- 
late z(5) for the inhomogeneous RTIM by considering the 
singular behavior of two quantities: First the probability 
distribution of the gap and then the asymptotic behavior 
of the surface autocorrelation function. 

From the distribution of the first gap one can de- 
duce the. dynamical exponent through the asymptotic 
relationfl 



lim 



In [P(lnei)] 
lnei 



1 



z(S) 



(4.16) 



As shown in Fig. |Uj the log-log plots of the distributions 
obtained for different decay parameters k, corresponding 
to relevant and marginal perturbations, have approxi- 
mately the same slopes. Consequently, the dynamical 
exponent z(8) is independent of the inhomogeneity. Sim- 
ilar results are obtained with the probability distributions 
of the second gap, which is expected to obey the same 
scaling relation as ig Eq. ( 4.16| ), however with an expo- 
nent z'{5) = z((5)/2.B As shown in Fig. Iljj, the probability 



Next we study the average surface autocorrelation 
function, the decay of which is connected to z(5) through 
Eq. (O). 
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FIG. 11. Average surface autocorrelation function in the 
Griffiths phase with the uniform distribution at ho = 3, 
A — 0.5 and for different values of k. The asymptotic slopes 
of the curves are given in parentheses. 



As shown in Fig. 
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the value of the decay exponent 
z(S) has only a weak dependence of the inhomogene- 
ity exponent k, both for relevant and marginal inhomo- 
geneities, due to the finite size effects which are expected 
to be larger when k is decreasing. Since the variation of 
z with the inhomogeneity is decreasing when the size of 
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the system increases, we expect that the surface autocor- 
relations are also asymptotycally independent of k. The 
value of z(S) is consistent with that obtained from the 
distribution of the gap and corresponds to the analytical 
bulk result in Eq.( |l.7| ). 

Thus we conclude that the presence of inhomogeneous 
disorder does not modify the form of the Griffiths-McCoy 
singularities and the dynamical exponent of the inhomo- 
geneous RTIM is the same as for the homogeneous one. 



V. DISCUSSION 

In this paper the influence of an inhomogeneous dis- 
tribution of the disorder in the vicinity of the surface of 
the RTIM has been studied. We considered a random 
version of the HvL model where, on the average, the per- 
turbation decays as a power l~ K of the distance I from 
the surface. The effect of the inhomogeneity was found 
to be relevant (irrelevant) for k < 1/2 (k > 1/2). 

The relevance-irrelevance criterion, which was shown 
to follow from a simple RW argument, can be generalized 
by making use of a more traditional approach Hj Near 
the bulk critical point, the charateristic length scale of 
the unperturbed system is set by the average correlation 
length £ ~ 5~" . The strength of the perturbation at a 
length scale I is measured by the mean change of the 

control parameter, 8 A (l) = l/ZEj=i Al ~ K ~ Al ~* '"■ This 
has to be compared, at the scale of the correlation length 
£, to the deviation from the critical point 6. Forming the 
ratio 5a(£,)/3 ~ A£~ K+1 ' V , we verify that the relative 
strength of the perturbation diverges (vanishes) at the 
critical point when k < \jv (k > 1/v). For the RTIM 
this is equivalent to the RW result of Section III since 
v = 2 in this case. 

Our next remark concerns correlated disorder in the 
RTIM, when the distribution of the couplings and/or 
transverse fields at different sites are algebraically cor- 
related with, for example, (Tri(J)ni +r (J)) ~ r - ". This 
type of correlated disorder has been recently studied in 
Ref. [2^. The relevance-irrelevance criterion has the same 
form as for inhomogeneous disorder with lo playing the 
role of k. Also the relation between relevant time and 
length s cale s, lnr r ~ C 1 "^, corresponds to the relation 
in Eq. (4.6) for relevant inhomogeneities. However the 
form of the critical singularities for the bulk and surface 
magnetizations are different in the two models. 

We close this Section with a discussion of a related 



problem, namely, the persistence probability, P pe r, for a 
RW in an inhomogeneous one-dimensional random en- 
vironment which is an inhomogeneous version of Sinai's 
model.u The transition probabilities, 7^ w i+i,i, 

are random variables and the corresponding probability 
distributions are inhomogeneous near the starting point 
of the walk, 1 = 1, deviating from the bulk distribu- 
tions by an aroftunt of order l~ K . We make use of an ex- 
act mapping^'Eil between the eigenvalue problem for the 



fermionic excitations of the RTIM in Eq. (2.4) and the 
eigenvalue problem for the master equation of the RW, 
with the following correspondences: 



Wl,l-i <-> hf 
m' 



-Fper(^) 
A 



\{L) 



(5.1) 



Here A m i n is the leading eigenvalue of the Fokker-Planck 
operator, L is the size of the system with free bound- 
ary conditions for the RTIM, whereas we take absorbing 
boundary conditions at both ends for the RW. Using the 
relations in Eq. (5.1), one can easily translate the results 



of Section IV. The inhomogeneity is relevant (irrelevant) 
for the diffusive and persistence properties of the RW for 
k < 1/2 (k > 1/2). In the case of a relevant inhomogene- 
ity with zero global bias, the persistence probability is 
either finite, for locally enhanced transition probabilities 
to the right, or it vanishes as 



[ p pcr(£)] av ~ exp (-const L 1 2h 



(5.2) 



for locally enhanced transition probabilities to the 
left. The relation between time and length scales, \nt, ~ 
i 1_K , leads to: 



[*per(*)] B 



exp 



-const (hit) 



(1-2k)/(1-k) 



(5.3) 



For marginal inhomogeneity the average persistence is 
given by 



[P pcr (i)] av ~(lui) 2e(A) 



(5.4) 



where the continuously vary ing exponent 9(A) is ob- 
tained as the solution of Eq. (|4.S|). 
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